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We study a model of a light-induced proton pump in artificial reaction centers. The model contains 
a molecular triad with four electron states (i.e., one donor state, two photosensitive group states, 
and one acceptor state) as well as a molecular shuttle having one electron and one proton-binding 
sites. The shuttle diffuses between the sides of the membrane and translocates protons energetically 
uphill: from the negative side to the positive side of the membrane, harnessing for this purpose the 
energy of the electron-charge-separation produced by light. Using methods of quantum transport 
theory we calculate the range of light intensity and transmembrane potentials that maximize both 
the light-induced proton current and the energy transduction efficiency. We also study the effect 
of temperature on proton pumping. The light-induced proton pump in our model gives a quantum 
yield of proton translocation of about 55%. Thus, our results explain previous experiments on these 
artificial photosynthetic reaction centers. 



I. INTRODUCTION 

It would be desirable to create an artificial system that 
exploits the basic principles of natural photosynthesis in 
order to produce energy in an usable form [TJ El [3j HJ O |6j 
7 . Indeed, natural photosynthetic structures efficiently 
convert the energy of light into chemical form [7J [H] . 

The overall energy transduction process in plant pho- 
tosynthesis occurs through a number of strongly coupled 
successive stages (see, e.g., [TJ EJ [5]). In the first step, 
light of the appropriate wavelength is absorbed by a light 
harvesting complex. The second step involves the conver- 
sion of electronic excitation energy to redox-potential in 
the form of the long-lived transmembrane charge separa- 
tion via multi-step electron transfer processes. The first 
two steps involve three constituents: (a) light-absorbing 
pigments, (b) an electron acceptor, and (c) an electron 
donor. In the third step, the energy stored in the electron 
subsystem is used for energetically uphill proton pump- 
ing, which generates the proton motive force across the 
membrane. 

The study of natural photosynthesis has inspired re- 
searchers to perform the photo-induced energy transduc- 
tion processes in the laboratory [TJ |5J |5[ |H [5] |u] [Tj $\ [Tfj] 
Hi] H2] . A convenient approach to photosynthesis in ar- 
tificial reaction centers is to use synthetic pigments, elec- 
tron acceptors and electron donors that are very similar 
in molecular structure to natural pigments (e.g., chloro- 
phylls, carotenoids and quinones). In this direction, the 
experimental model proposed in Refs. [HI HHj provides 
a paradigm for the conversion of light energy to a pro- 
ton potential gradient. These seminal works 0110] have 
motivated research in the design and synthesis of new 
artificial photosynthetic systems [T2J [21 [H] (i.e., light- 
harvesting antennas and reaction centers) and triggered 
considerable experimental |T6j [TTl [TBI ESI EH] and theoret- 
ical EJJ [22j E3J EI] activities to investigate more sophis- 
ticated and more efficient mechanisms for the conversion 
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FIG. 1: (Color online) The top figure presents the triad 
(donor "D" , photo-sensitive part "B,C", and acceptor "A") 
and the shuttle "S" [2] ITO] . These are enclosed by color cir- 
cles, which are schematically shown in the bottom figure. 
The tetraarylporphyrin group acts as a photosensitive moi- 
ety (B,C) (inside the green circle in the top structure). This 
is connected to both a naphthoquinone moiety fused to a nor- 
bornene system with a carboxylic acid group (which acts as 
an electron acceptor (A)) and to a carotenoid polyene (which 
acts as an electron donor (D)). 2,5-diphenylbenzoquinone is 
the proton shuttle (S), denoted by a pink hollow circle in the 
structure and by a solid pink circle in the cartoon. 



of light energy. 

The transformation of light energy into the electro- 
chemical gradient of protons across the membrane can 
be quantitatively characterized by the quantum yield (or 
quantum efficiency), of proton translocation. This 
parameter is defined as the total number of translocated 
protons divided by the number of photons absorbed by 
the triad [9]. A quantum yield of the order of 0.4% 
has been measured in Ref. [9]. A much higher quan- 
tum efficiency, ~ 7%, for the conversion of photons 
into ATP molecules, was found in Ref. [10]. As argued 
in Ref. [TO] , the actual quantum yield of ATP formation 
could be of the order of 15%, if we take into account the 
real rate of light absorbance, which is ~ 50%. Near four 
protons are necessary for the synthesis of a single ATP 
molecule. This means that the real quantum yield $ of 
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FIG. 2: (Color online) Schematic diagram of the light-induced 
proton pump across the lipid bilayer in a liposomic membrane. 
A molecular triad D-BC-C is symmetrically inserted in the 
lipid bilayer. The different stages in the proton pumping pro- 
cess are here denoted by (a,b,c,d,e,f). The two bluish vertical 
rectangles on both sides schematically represent two proton 
reservoirs with electrochemical potentials fip and /un- These 
two proton reservoirs correspond to the aqueous phases in- 
side and outside of the liposome, respectively. The shuttle 
molecule S, is shown as a pink-colored oval and the proto- 
nated neutral shuttle is shown as a yellow oval. This shuttle 
freely diffuses in (d) (the black scribbled curves represent the 
thermal stochastic motion of the shuttle) across the mem- 
brane to transport a proton from the lower proton potential 
/in to the higher proton potential /jp side of the membrane, 
where (fip—fj,^) denotes the total potential difference between 
the two reservoirs. 



proton translocation measured in Ref. [TU] can be about 
60%. The total thermodynamic (or power-conversion) 
efficiency, 77, of the light-to- ATP conversion process is 
estimated in Ref. [TO] as r] ~ 4%. 

In the present paper, using methods from quantum 
transport theory [251 ED 123 HE] , we analyze the photoin- 
duced electron and proton transfer in a molecular triad 
inserted into a liposomal membrane, which contains a sin- 
gle molecular shuttle. We calculate the photon-to-proton 
quantum yield $ ~ 55% (and the thermodynamic effi- 
ciency 77 ~ 6.3%) for the resonant tunneling conditions, 
when the reorganization energy, A, of the electron transi- 
tions matches the detuning 5 between the electron energy 
levels: A ~ 5. 

We note that due to a small optimal value of the re- 
organization energy (A ~ 400 meV) the charged recom- 
bination process in the triad is described by the inverted 
region of the Marcus formula [20J [30] El • This further 
enhances the performance of the system. Our results 



explain experiments made in Ref. [10] using artificial 
photosynthetic centers. The obtained power-conversion 
efficiency corresponds to the highest value, r\ ~ 6.5%, 
achieved recently with polymer solar cells [35]. It is ex- 
pected that the proton current and the efficiency should 
increase with increasing the number of the shuttles in the 
membrane. 

This article is organized as follows. In Sec. II (see 
also the Appendix) we introduce the basis set for the 
system and write the Hamiltonian of the problem. In 
Sec. Ill, we present the master equation for the density 
matrix coupled to the Langevin equation describing the 
diffusive motion of the shuttle in the lipid bilayer. In Sec. 
IV, we numerically solve these equations and analyze the 
light-induced proton pumping process. In Sec. V we 
summarize our results. 



II. MODEL 

We use a slightly modified version of the well-accepted 
model already presented, e.g., in Refs. [TU]. In this 
model the reaction center is a molecular triad containing 
an electron donor and an electron acceptor both linked 
to a photosensitive porphyrin group (shown in Fig. 1). 
The triad molecule (D-BC-A) is inside the bilayer of a 
liposome. The lipid bilayer also contains freely diffus- 
ing 2,5 diphenylbenzoquinones, acting as proton shut- 
tles. The molecular triad absorbing a photon establishes 
a negative charge near the outer surface and a positive 
charge near the inner surface of the liposome, by gener- 
ating charge separated species D+-BC-A - . The freely 
diffusing quinone shuttle translocates an electron-proton 
pair across the membrane and neutralizes the molecular 
triads. In Fig. 2 we schematically illustrate the process 
of light-induced proton pumping in liposomes by artificial 
photosynthetic reaction centers [9j [10] • The transmem- 
brane proton pumping requires a symmetric arrangement 
of the molecular triad (of length ~ 8 nm) inside the bi- 
layer and with a specific direction: with the acceptor (A) 
site towards the outer membrane of the liposome (the 
negative (N) side of the membrane) , and with the donor 
(D) towards the inside of the liposome (the positive (P) 
side of the membrane) [9] [10] ■ 

The energy diagrams of the electron and proton sites 
are shown in Figs. 3 and 4. There are two electrons in the 
system, one of which is initially on the D site, and an- 
other electron is on the lower energy level B. The quinone 
molecular shuttle has one electron state S (denoted by S, 
instead of S G ), and one proton state Q (denoted here by 
Q instead of S p ). Thus, S denotes the shuttle electron 
state and Q denotes the shuttle proton state. 

The overall process leading to the proton translocation 
from the N-reservoir with a lower proton potential, /jn, to 
the P-reservoir with a higher electrochemical potential, 
fip , can be considered as a sequence of eight stages (most 
of which are shown in Fig. 2). 

• Step I: The photosensitive moiety of the molecular 
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3: (Color online) Energy diagram depicting the en- 
ergy levels of states involved in an artificial photosynthetic 
reaction center, before the diffusion of the shuttle to the P- 
reservoir. The subfigures (a,b,c) correspond to the stages 
(a,b,c) in Fig. 2. The left and right panels represent electron 
and proton energy levels, respectively. The abbreviations D, 
B, C, A, S are the same as used in the text and in Fig. 1. 
Also, and xa represent the spatial coordinates of the sites 
D and A, respectively. The thick brown arrows denote the 
path the electrons follow in this energy diagram, generating 
charge separation, in (b), and shuttle charging and protona- 
tion in (c). Initially, light excites an electron from B to C, 
and eventually to A, making it A - . Afterwards, in (b), the 
donor D loses an electron, thus becoming D + , and that elec- 
tron moves to BC. Later on, the shuttle S in (c) receives the 
electron from A. 



triad absorbs light and an electron goes from the 
ground state B to the excited state C (see Fig. 3b). 

• Step II: The unstable excited state C transfers the 
electron to the acceptor A, producing an unstable 
charge-separated intermediate species D-BC + -A~. 

• Step III: The unstable intermediate charge- 
separated species is rapidly rearranged to a rela- 
tively stable charge-separated form (D+-BC-A - ) 
by the thermal electron transfer from the state D 
to the state B + having a lower energy than the state 
D (Fig. 3b). 

• Step IV: The shuttle in the position near the N-side 
of the membrane accepts an electron from A~ and 
becomes negatively charged. 

• Step V: The shuttle molecule receives a proton from 
the N-reservoir and becomes neutralized (Fig. 3c 
right panel). 
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FIG. 4: (Color online) Energy levels involved in an artifi- 
cial photosynthetic reaction center. This figure is similar to 
Fig. 3, but now the energy profile corresponds to the stage 
after the shuttle diffuses to the P-reservoir. Here the subfig- 
ures (d,e,f) correspond to the stages (d,e,f) in Fig. 2. The left 
and right panels represent proton and electron energy levels, 
respectively. The thick brown arrows denote the path fol- 
lowed by the electron (e) and proton (f). In (d), an electron 
on the shuttle S moves to the donor site D, neutralizing it in 
(e). This electron transition in the right panels increase the 
proton energy of the shuttle, as shown in the left panels (from 
(d) to (e)). The proton finally leaves the shuttle in the left 
panel of (f). 



• Step VI: The neutral shuttle slowly diffuses through 
the lipid bilayer and carries the electron and the 
proton to the P-side of the membrane and to the 
D-site (stage (d) in Fig. 2). 

• Step VII: The shuttle gives away the electron to the 
positively charged site D + (stage (e) in Fig. 2 and 
Fig. 4e). 

• Step VIII: The shuttle is dcprotonated by donating 
the proton to the P-reservoir (Fig. 4f). 

This sequence of eight steps describes the photo-induced 
electron transfer that generates the intra-membrane re- 
dox potential, which in turn drives the energetically up- 
hill vectorial translocation of protons by the shuttle. 

Electrons in the states i (= D,B,C,A,S) and protons in 
the state Q are characterized by the corresponding Fermi 
operators af , a, and 6q , &q , with the electron population 
operator rii and the proton population hq. We assume 
that each electron or proton state can be occupied by a 
single electron or a single proton. Spin degrees of freedom 
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are neglected. The proton site on the shuttle, denoted by 
Q, can be populated from the N-reservoir provided the 
shuttle is within the transition length Lq from the re- 
side of the membrane. The protonated shuttle, located 
within the transition (or tunneling) range from the P- 
side of the membrane, can donate its proton to the P- 
reservoir. Protons in the reservoirs are described by the 
Fermi operators d^ a ,dk, a , where a = N,P; and k is an 
additional parameter which has the meaning of a wave 
vector in condensed matter physics [25, 26, 27, 28]. The 
number of protons in the reservoirs is determined by the 
operator J2k N ka, with N ka = d^ a d ka . 



A. Hamiltonian 

The Hamiltonian of the electron-proton system, 

H = H + H dil + H tI + H B , (2.1) 

has a term H related to the energies E^ of the electron 
eigenstates (i = D,B,C,A,S), and to the energy eq of a 
proton, on the shuttle: 

H = EjUj + eQUQ + Mdb(1 - np)(l - n B - Tip) 

i 

- u d a(1 - n B )n A - uba(1 ~ n B - «c)«A 



usq n s uq. 



(2.2) 



We include here the electrostatic interaction between the 
electron sites, «db!«da,"ba, and the Coulomb attrac- 
tion usq between the electron and proton sites on the 
shuttle. It is assumed that the empty donor state D 
(with rip = 0) as well as the empty photosensitive group 
B and C (ne + nc = 0) have positive charges, and 

«DB = "DC, "CA = «BA- 

The term, 



Ha 



-A DB a D a B - A A c a A a c - A DS (a;) a D a s 



- A A s (x) a A a s - F(t) a B a c + h.c. 



(2.3) 



describes the tunneling of electrons between the sites D- 
B, C-A, A-S, and D-S, with the corresponding ampli- 
tudes Ajj/. Notice that the tunneling elements Ads (a;) 
and Aas(^) depend on the shuttle position x. The 
Hamiltonian is also responsible for the electron tran- 
sitions between the states B and C induced by the elec- 
tromagnetic field (light), F(t) — F exp(iw i), with a 
frequency uq and an amplitude Fq. Proton transitions 
between the shuttle (site Q) and the N- and P-proton 
reservoirs are governed by the Hamiltonian 

fftr = - E T ^ x ) d L b Q ~ E T *^ x ) b Q d ^ (2-4) 



with the position-dependent coefficients, Tk a (x). We 
have chosen the following form of Tk a (x): 

TkN(x) = T kN 6[x - (x N - Lq)], 
Tkp(x) = T kP 0[x P + L Q - x], 



where 0(x) is the Heaviside step function, and the param- 
eter Lq defines the proton loading range of the shuttle. 

B. Interaction with the environment 

To take into consideration the effect of a dissipative en- 
vironment we consider the well-known system-reservoir 
model [251 133 133], where the medium surrounding the 
active sites is represented by a system of harmonic oscil- 
lators with the Hamiltonian: 



E 



(2.5) 



where Xj , pj are the positions and momenta of the oscil- 
lators with effective masses rrij and frequencies ujj. The 
parameters Xu determine the strengths of the coupling 
between the electron subsystem and the environment. 
The system of independent oscillators are conveniently 
characterized by the spectral functions Jhi(lo), defined 

by 

, hAu) = £ rn 3 ^( Xjl -x 3l ,f ^ _ ^ (2 g) 



so that the reorganization energy A^< , related to the i 
i' transition, has the form 



Xw 



Jw \v) — 2 ■ ' > 

j 



UJ 



With the unitary transformation U — Y[ t Ui, where 



Ui = exp 



y^PjXjjnj 
j 



(2.8) 



we can transform the Hamiltonian H to the form H' = 
HU, becoming (after dropping the prime) 

h = ff -E A «' e(i/2)(ei " ?D 4^ 

ii' 

- F(t)e^™ B ~^ a B a c - F*(t) a\. a B e^^-Sc) 



k(x 



ka 



y - 



2 2 

mjLuljXj 



2m, 



(2.9) 



where a = N,P, and the tunneling coefficients, A*^ = 
Aj/j, take non-zero values only for transitions between 
the sites D and B, A and C, A and S, as well as D and 
S. The stochastic phase operator & is given by 



(2.10) 
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The result of this transformation follows from the fact 
that, for an arbitrary function <fr(xj), the operator U 
produces a shift of the oscillator positions: 



( l x- \ 

I x j ' 2 2-^i x i i ' rt " i I 



(2.11) 



This transformation also results in the phase factors for 
the electron amplitudes (see Eq. (9)). 

The basis sets, composed of the electron-proton eigen- 
states, and their corresponding energy eigenvalues are 
presented in an Appendix. Thus, the reader is encour- 
aged to read this short Appendix before proceeding fur- 
ther. 



III. TIME EVOLUTION OF DENSITY MATRIX 

A. Master equations 

To describe the time evolution of the diagonal elements 
of the density matrix, (p m ), we write the Heisenberg 
equation for the operators p m with the subsequent av- 
eraging over the environment fluctuations and over the 
states of the proton reservoirs: 



(p m ) = -(i[p m ,H d i T ]-) - (i[p m ,H iT ]-). 



(3.1) 



The protons in the reservoirs (a = N,P) are character- 
ized by the Fermi distributions, 



F a (Eka) 



exp 



/ Eka — Pa 



T 



(3.2) 



with the temperature T (fee = 1). The electrochemical 
potentials /in and pp, correspond to the negative (N) and 
positive (P) proton reservoirs, respectively. The proton 
motive force (Ap) across the membrane is given by 



Ap = p P - p N = V 



2.3 RT 
F 



(ApH) 



(3.3) 



where R and F are the gas constant and Faraday con- 
stant, respectively, and V is the transmembrane voltage 
gradient. Hereafter we change Ap by changing the pH 
of the solution by ApH. 

The contribution of the transitions between the shuttle 
and the proton reservoirs to the time evolution of the 
density matrix is described by the second term in the 
right hand side of Eq. (12), which can be calculated with 
methods of quantum transport theory |25l 126) 

(i[ Pm ,tf tr ]_) = J2 [7nm(x){Pm) - 7™„(z)(p»>] , (3-4) 
n 

with the relaxation matrix 



7 tr 
/ mn 



(x) = y^r a (x) {|&q, m n| 2 [l - F a (w nm )] 



\ 2 F a (uj mn )}. (3.5) 



Here we introduce the frequency-independent coeffi- 
cients, 

T a {x) = 2tt£ \T ka (x)\ 2 S(u - E ka ), (3.6) 

k 

which determine the transition rates between the shut- 
tle state Q and the sides of the membrane (N- and P- 
reservoirs). Notice that these coefficients are functions 
of the shuttle position x. 

The transitions between the electron levels are de- 
scribed by the Hamiltonian iJdir, which can be written 
as 

Hdir ^ Aran Pm,n ^ Pn.m ^mni (3-7) 

ran mn 

with the functions 

A mn — QDB(aB a o)mn + Qc\{(^x a c)mn + QsA(aA a s)m« 

+ QsB(al,as)mn + QcBia^a^mn, (3.8) 

which arc defined as superpositions of the heat-bath op- 
erators 

Qu> = A</<exp[(*72)(6 -iv)\ 

= A i > i e3qp[(i/2)^2p j (t)(x ji -x ji f)], (3.9) 

3 

for the pairs of the electron sites (ii') = 
(DB),(CA),(SA),(SD), whereas for the pair (CB) 
we have 

Qcb = F exp(iuj t) exp[(i/2)^2p j (t)(x jC ~ x jB )], 

(3.10) 

In the case of a high-enough temperature of the bath [30] , 
the cumulant functions of the unperturbed operators QW 
are determined by the relations: 

(3.11) 

The contribution of the electron transitions to Eq. (12) 
is determined by the term 

( ^[pWij -^dir] — ) ^ ^ {•AmnPmn ~ •A-nmPnm) H~ h.C. 

n 

(3.12) 

Within the theory of open quantum systems developed in 
Refs. [28 , the correlation function (A mn p mn ) is propor- 
tional to the density matrix elements of the system, (p m ), 
with coefficients defined by the unperturbed correlators 
(Amn(t),Amn (*')) °f the bath operators: 



(A mn (t)pmn(t)) = i J dheit-hy"™^-^ 

x {{AM(t),A<M(ti)){p m (t)) 
- (A^it^A^itMp^t))}, 



(3.13) 
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where 



<Q^(*),Q^ t (*l)>l(4oc)mn| a 
(Qdb(*)' ( 3DB t ( i l))K a B a D)m™| 2 
(QcAW^CA(<l)>K«Vc) m J 2 

(4°DW J QsD t (*l)>K«D«s) mn | 2 , 

(3.14) 



and the reverse expression can be obtained for the cor- 
relator (.A^n (*i),^mn(*)>- Thc formula (24) is valid in 
the case of weak tunneling and weak driving force Fq. 
The effects of quantum coherence are also neglected here. 

Finally, we derive the master equation for the density 
matrix of the system, 

(x)(p n ), (3.15) 

n n 

with the total relaxation matrix 

+ (fSA)mn + (KSD)ran + (^CB)mn, (3.16) 

containing the contribution of proton transitions to and 
from the shuttle, 7^(1), together with the Marcus rate 
(kcbJuih describing the light-induced electron transfer 
between the sites B and C: 



(«Bc)r 



\F a 



x exp 



\{a B a c ) mn \ 



{u mn + ujq + XbcY 
4A BC T 



x exp 



(Umn — WQ + ^BcY 

4A BC T 



(3.17) 



as well as the rates related to the electron transfers be- 
tween the pairs of sites (ii') = (DB),(CA),(AS), and 
(DS): 



{j^ii' )mn — | ^i'i \ 



\(a\,ai) mn \ 2 + \{a\,ai) nm Y 



x exp 



1A„V 



(3.18) 



We note that the tunneling coefficients Aas and Ads 
depend on the shuttle position x. 



by the overdamped Langevin equation for the shuttle po- 
sition x: 



dx 
~dt 



dU{x) 
dx 



(3.19) 



Here ?7drag is the drag coefficient of the shuttle in the lipid 
membrane, and the thermal fluctuation of the medium is 
modelled by a zero-mean delta-correlated Gaussian fluc- 
tuation force C(t), (C(t)) = 0, 



(C(*)C(0>=2Wra(*-*'), 



(3.20) 



where T is the temperature of the medium (fes=l). 
The diffusion of the shuttle is determined by the diffu- 
sion coefficient D s = T/?7drag- The potential U{x) in 



Eq. (3.191 is responsible for the spatial confinement of 



the hydrophobic shuttle (quinone) inside the lipid mem- 
brane. 



IV. RESULTS AND DISCUSSIONS 

To analyze the light-induced proton pumping process 
quantitatively, we use the standard Heun's algorithm to 
numerically solve the twenty coupled master equations 
(26) along with the equation (30) for the shuttle. For ini- 
tial conditions we have assumed that at t = 0, pi t ± = 1, 
and the other elements of the density matrix are zero 
(this corresponds to one electron on site D and another 
electron on site B with no electrons and no protons on 
the shuttle). We also assume that at t = thc shuttle 
is located nearby the acceptor (A): x(t = 0) = xa — %n- 
Throughout our simulation we focus on the long-term 
asymptotic regime, where the effects due to the influ- 
ence of transient processes have been smoothed out. The 
time-homogeneous statistical properties are obtained in 
the long-time limit after the temporal and ensemble av- 
eraging are performed. 

Thc efficiency (quantum yield) of the proton pumping 
device is defined by the formula: 



number of protons pumped 
number of photons absorbed 



The photon absorption rate, Kb— >d is approximately 
equal to the rate of light-induced transitions from the 
state B to the state C. Thus we assume, 



KB^C 



(4.1) 



where I p is the proton current (the number of protons, 
A p , translocated across the membrane per unit of time). 



B. Equation of motion for the shuttle 



Diffusive motion of the shuttle in the lipid 
bilayer 



We assume that the shuttle moves along the linear 
molecular triad (Fig. 1), and this motion can be described 



In Fig. 5 we present the diffusive motion (see Fig. 5(a)) 
of the shuttle in the lipid bilayer together with the time 
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FIG. 5: (Color online) (a) Stochastic motion of the shuttle 
with time. The horizontal black dashed lines denote the bor- 
ders of the membrane, x-^ = 40 A, xp = - 40 A. Via this 
diffusion the shuttle transports protons and electrons through 
the membrane, (b) Variation of the electron and proton pop- 
ulation on the shuttle. Note that the proton density (red 
curve) and the electron density (black curve) mostly coincide 
in (b). (c) Number of protons pumped versus time. The 
main parameters used here are the light intensity I = 0.138 
mWcm -2 , temperature T — 298 K, and the chemical poten- 
tials /j,p = 110 meV and /^n = —110 meV. The light intensity 
I corresponds to the photosensitive BC-group with a dipole 
moment ~ |e| x 1 nm, where e is the electron charge. 



energy) of the order of 0.3 eV [51] . Even smaller reor- 
ganization energies (A ~ 230 meV) have been measured 
for the porphyrin-fullcrene dyads [35 . The initial stages 
of electron transfer in bacterial reaction centers are 
also characterized by a low reorganization energy: A ~ 
70-300 meV, depending on the environment. This is 
due to the fact that the bacteriochlorophyll molecules 
(and the molecules of porphyrin involved in our molecu- 
lar triad) contain highly delocalized 7r-electron systems. 
In the next section, we also analyze how sensitive the 
results are to changes in the values of A. 

Electrochemical measurements [1J show that the 
energy of the carotene(D)-porphyrin(BC)-quinone(A) 
molecular triad sweeps from the value ~1.9 eV (the 
first excited state of the porphyrin, D-B 1 C-A), to the 
energy ~1.4 eV, related to the intermediate state D- 
BC + -A~, and, finally, to the energy, ~1.1 eV, of the 
charge-separated state D+-BC-A". We assume here 
that the energy of the first excited state of the porphyrin, 
Ec — Eb, is 1908 meV, which corresponds to a photon 
wavelength of 650 nm as used in experiments [9] [10]. We 
have taken the energy gap between the site C and A 
to be approximately equal to the reorganization energy, 
(Eq ~ Ea) ~ A = 400 meV. This gap is about the en- 
ergy difference between the D-B 1 C-A and D-BC + -A~ 
states. 

The energies of the electron sites S and A are compa- 
rable, (Ea — Eg) — 300 meV, due to a structural simi- 
larity of the quinone shuttle (S) and quinone moiety of 
the molecular triad. The protonation of the shuttle leads 
to the lowering of the electron energy on site S due to 
the electron-proton Coulomb attraction [9], usq ~ 360 
meV. The other Coulomb interaction terms are chosen as 
«db = uba = 120 meV and «da = 60 meV. These values 
correspond to the electrostatic interaction of two charges 
located at distances 4 nm and 8 nm, respectively (in a 
medium with a dielectric constant ~ 3). Furthermore, 
we assume that Ed — Eq = 400 meV and £q = 200 meV. 
We have chosen eQ such that, for the above mentioned 
parameters, the device works well at the transmembrane 
potential difference ~ 200 mV. 

We choose fip = 110 meV, = —110 meV, the res- 
onant tunneling rates A/H = 15 ns , T/h = 1.5 ns _1 , 
and the reorganization energy for the light-induced elec- 
tron transfer, Abc ~ 80 meV. The majority of parameters 
in our model are deduced from experimental data. The 
rates of electron transfer reactions are given by 



dependencies of the electron and proton populations of 
the shuttle (Fig. 5(b)) complemented by the time evolu- 
tion of the number of pumped protons (Fig. 5(c)). We 
assume that the reorganization energies for the thermal 
electron transfers are low enough to provide a high per- 
formance of the system: 

A ~ A DB ~ A A c ~ A A s ~ A DS ~ 400 meV. 

This value of A is quite common for porphyrin-quinonc 
dyads having a lower limit (the internal reorganization 



A - K D ^ B ~ 26 (IS , 
KA^S — KS^A ^ 20 ^S" 1 . 

Therefore, the loading and unloading time scales of the 
shuttle are about 0.05 fis. The shuttle has enough time 
to be loaded and unloaded with electrons and protons 
when it enters the loading/unloading domain with a size 
about the electron tunneling length, Ltun ^0.5 nm, and 
the proton transition length, Lq ~0.2 nm. Figures 5(a,b) 
show a time synchronization between the spatial motion 



8 



of the shuttle and the time variations of the shuttle pop- 
ulations. 

It follows from Fig. 5(c) that in 1 ms the shuttle per- 
forms near 16 trips and translocates 10 protons through 
the membrane, provided that the light intensity I — 0.133 
mWcm~ 2 . We assume that the diffusion coefficient D s is 
of the order of 2 nm^s" 1 [37 , and the dipole moment 
of the BC moiety is about |e| xl nm, where e is the elec- 
tron charge. The number of photons absorbed in 1 ms 
is ~ 18. Thus, the approximate quantum yield of the 
pumping process is ~ 55 %. For this parameters, the dif- 
fusive motion of the shuttle is the slow and rate-limiting 
step of the pumping process. 

B. Robustness of the model 

To show a tolerance of the system to variations of 
parameters we explore here the parameter space of our 
model. Keeping fixed the energy difference between the 
sites B and C, we calculate and plot (see Fig. 6) the 
pumping efficiency $ (photon-to-proton quantum yield) 
as a function of the reorganization energy, A, and the 
energy gap, 8, 

A ~ Adb ~ Aac ~ ^as ~ <^DS, 
8 = Ec — Ea = Eg — -Ed , 

between the energy levels Ec and Ea, and between the 
levels Eg and E&. The figures 6(a), 6(b), 6(c) corre- 
spond to the different values of detuning between the 
acceptor energy level, Ea, and the electron energy level 
on the shuttle, E s : E A - E s = 100 meV (Fig. 6(a)); 
Ea- E s = 300 meV (Fig. 6(b)); E A - E s = 500 meV 
(Fig. 6(c)). These plots clearly demonstrate the exis- 
tence of quite wide areas in the plane A - 8, where the 
pump performs with maximum efficiency. For the detun- 
ing Ea—E$ = 100 meV (Fig. 6(a)) the pumping efficiency 
reaches its maximum, $ ~ 48%, in the region of param- 
eters (in meV): 270 < A < 500, and 400 < 5 < 700. In 
this region, the energy gaps between the redox sites are 
close to the reorganization energy, which results in higher 
site-to-site tunneling rates and, consequently, in a high 
pumping efficiency. 

The higher pumping efficiency, $ ~ 55%, can be 
achieved at the detuning Ea — E$ = 300 meV (Fig. 6(b)). 
In this case the parameter 8 can be tuned in such a way 
that the energy gaps between all relevant electron sites 
are equal to the reorganization energy: 

(E c — E A ) ~ {Ea - E s ) 

~ (E s ~ u S q - Ed) ~ (E D - E B ) ~ A. 

We recall that a shuttle populated with a proton has 
the electron energy, Es — u$q, which differs from the 
initial value Eg by the charging energy, wgQ ~ 360 meV. 
Summing all the above-mentioned detunings and taking 
into account the energy difference, Eq — Eb — 1908 meV, 
between the optically-active levels B and C, we estimate 




FIG. 6: (Color online) Contour plots presenting the variations 
of the quantum efficiency <& with the reorganization energy A 
and with the energy gap 5, where 8 = Ec — Ea = Es — 
Eo- The parameters used here are: light intensity / = 0.138 
mWcm~ 2 , temperature T — 298 K, and chemical potentials 
fip = 110 meV and /in = —110 meV. The detunings take the 
following values: (a) E A ~E S = 100 meV; (b) E A -E S = 300 
meV; (c) E A - E s = 500 meV. 



the optimum values of the reorganization energy A and 
the detuning 5: 



A ~ 5 ~ (E c - E B - u S q)/4 = 387 meV. 



The maximum of the efficiency in Fig. 6(b) is observed 
at S ~ A ~ 400 meV, which is very close to our estima- 
tions. For a larger energy gap, Ea — Es — 500 meV (see 
Fig. 6(c)), the proton pumping efficiency $ decreases and 
the region of the optimum parameters shrinks compared 
to Fig. 6(b). 
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C. Effects of the resonant tunneling rates 

The fine-tuning of tunneling couplings between active 
electron sites is feasible in some nanostructures. This 
tuning can be implemented by changing the site-to-site 
distance, as well as by varying the height of the potential 
barriers (see, e.g., [38]). Artificial photosynthetic sys- 
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FIG. 7: Color online) Proton pumping quantum efficiency 
<3> versus resonant tunneling rate A, at different reorganiza- 
tion energies A, shown in (a,b), and for different detunings S, 
shown in (c,d). Note that A here represents A/ft, since we set 
ft = 1. We use the following parameters: / = 0.138 mWcm -2 , 
T = 298 K, hp = 110 meV, u. N = -110 meV, and the energy 
gap (E\ — Es) — 300 meV. Panels (a,b) are plotted at fixed 
S — 400 meV, whereas in (c,d) the reorganization energy is 
fixed, with A = 400 meV. 

terns, such as the molecular triads, also allow to engi- 
neer desirable tunneling and electrostatic properties of 
the structures [T] with the goal to achieve the highest 
possible efficiency. As in colloidal nanocrystals [38], this 
can be done by inserting additional molecular bridges be- 
tween the side centers D, A and the photosensitive part 
BC utilizing the exponential dependence of electron tun- 



neling rates on the distance [S]. 

In Fig. 7 we illustrate the variation of the proton pump- 
ing efficiency $ as a function of the resonant tunneling 
rate (A/7i) for different values of the reorganization en- 
ergy A (Figs. 7(a), 7(b)) and the energy gap 6 (Figs. 7(c), 
7(d)). The detuning, Ea — Eg, is fixed to the value 300 
meV for all plots in Fig. 7. 

In Fig. 7(a) we plot four curves, ^(A), for the following 
set of reorganization energies: A = 100, 130, 200, 400 
meV and for a detuning 5 = 400 meV. In Fig. 7(b) the 
efficiencies ^(A) are plotted for the reorganization en- 
ergies: A = 500, 800, 1000, 1200 meV, for the same 
detuning S. Similar dependencies, $(A), are depicted in 
Fig. 7(c) for S = 100, 130, 200, 400 meV and in Fig. 7(d) 
for S = 500, 600, 700, 800 meV. In both, Figs. 7(c) and 
7(d), the reorganization energy A is equal to 400 meV. 

It follows from Fig. 7 that initially the proton pumping 
efficiency rapidly increases with increasing A, followed by 
its saturation for higher values of the resonant tunneling 
rate. The saturation limit depends on the reorganization 
energy A as well as on the energy gap <5. For the optimum 
values of A and 5: A ~ 6 ~ 400 meV, the pumping effi- 
ciency is sufficiently high, $ ~ 55 %, even for moderate 
tunneling rates, A/fi < 5 ns _1 . 



D. Effects of Coulomb interactions 

In Fig. 8 we plot the efficiency 4> versus the dielec- 
tric constant e of the medium, to explore the effects of 
the Coulomb couplings udb, "BA, and «da on the perfor- 
mance of the proton pump. The electrostatic interactions 
between the photosensitive part B and C and the donor, 
mdb, between the sites B and C and the acceptor, mba, 
and between the donor and the acceptor, it da, are in- 
versely proportional to the dielectric constant e and to 
the distance between the relevant sites. For example, we 
have 

e 2 

UDB 



47r e £ r DB 



where tdb characterizes the spatial separation of the 
sites D and B, and £o is the vacuum permittivity. The 
Coulomb interactions between the sites D and B and be- 
tween the sites B and A (with tdb = ^ba = 4 nm) are 
decreased from 360 meV to 36 meV when the dielectric 
constant e scans the range from 1 to 10 (see Fig. 8(a)). 
We note that in our model rrjA — 8 nm, so that «da — 
mdb/2. Figure 8(b) shows the efficiencies for differ- 
ent values of the detuning 6:6 = 100, 200, 400, 600 
meV, for A = 400 meV. Moreover, in Fig. 8(c) we plot 
the efficiencies, <&(£), for fixed detuning 6 = 400 meV and 
for A = 100, 200, 400, 600 meV. 

It should be emphasized that near the optimum work- 
ing point (at A ~ 6 ~ 400 meV) the pump operates with 
the high efficiency, $ ~ 55 %, which practically does not 
depend on the dielectric properties of the medium. 
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FIG. 8: (Color online) (a) Coulomb energies mda and mdb 
versus the dielectric constant e of the medium, (b) Proton 
pumping efficiency $ versus dielectric constant e for different 
values of S and for A = 400 meV. (c) The pumping efficiency 
$ as a function of the dielectric constant e for different reor- 
ganization energies A and at the fixed detuning 6 = 400 meV. 
The other parameters are the same as in Fig. 7: I = 0.138 
mWcm" 2 , T = 298 K, /i P = 110 meV, p N = -HO meV, and 
E A - E s = 300 meV. 



E. Effect of light intensity 

In Fig. 9 we plot the proton current as a function of the 
light intensity for different values of the temperature. At 
zero light intensity the proton current is zero. Initially, 
with increasing light intensity, the proton current also in- 
creases linearly and then saturates around 0.2 rnW cm -2 . 
This saturation is probably caused by the slow diffusion 
of the shuttle inside the lipid membrane. A similar in- 
tensity of saturation (~ 0.1 mW/cm 2 ) has been observed 
in experiments 10 . 

In a warm environment, the shuttle moves faster and 
carries more protons. To do this, the system should ab- 
sorb more photons, so that at high temperatures a full 
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FIG. 9: (Color online) (a) Proton current versus light inten- 
sity I for different temperatures, at /un = —110 meV, and 
fip = 110 meV. Notice that the proton current is roughly lin- 
ear for small intensities of light, but it saturates with higher 
light intensity. In this saturation region, the proton current 
is larger with higher temperatures, (c) The standard devia- 
tion, cjp, of the number N p of pumped protons as a function 
of the light intensity I, for different temperatures, (b) The 
pumping quantum efficiency $ decreases with light intensity 
for all temperatures shown. 



saturation takes place at higher light intensities. We note 
that the low saturation limit obtained above and mea- 
sured in the experiment |10j with the carotene-porphyrin- 
quinone triads is far below the average intensity of solar 
light, I ~ 30 mW/cm 2 [B]. This fact points to the relative 
inefficiency of the energy-conversion process available at 
normal daylight conditions. An ideal highly-efficient pho- 
tosynthetic system should not have any saturation limits 
for the standard daylight intensity of light. 



11 



It is evident from Fig. 5 that the number of protons, 
N p , translocated across the membrane fluctuates in time. 
To estimate these fluctuations we calculate the standard 
deviation, 
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= yJ(N*) (N p )*, 



(4.2) 



which characterizes the magnitude of its shot noise. The 
dependence of the noise level cr p on the light intensity / 
is shown in Fig. 9(b). For an intensity of light / ~ 0.14 
mW 2 /cm 2 (when ~ 10 protons are translocated across 
the membrane and the efficiency, $ ~ 55%, is sufficiently 
high) the uncertainty a p in the number of pumped pro- 
tons is about 1.3. In Fig. 9(c) we demonstrate that the ef- 
ficiency $ of the light-induced pumping decreases mono- 
tonically with increasing light intensity. At low light in- 
tensities, a relatively small number of photons are ab- 
sorbed per unit time. Thus, a higher fraction of the 
absorbed photons is used for the uphill pumping of the 
protons. 



F. Effect of temperature 

Figure 10 shows the effects of temperature on the 
pumping current and on the efficiency of the photosyn- 
thetic device for different values of the light intensity. 
The temperature effects appear in the light-induced pro- 
ton pumping dynamics through two factors: (i) The elec- 
tron transfer rates, including the loading and unloading 
rates of the shuttle, increase with increasing temperature, 
(ii) The diffusion coefficient of the shuttle increases with 
temperature. Because of this, the shuttle can perform a 
higher number of trips to translocate protons at higher 
temperatures. Here the electron transfer reactions are 
not rate-limiting ones. The diffusive trips of the shut- 
tle from the N terminal to the P terminal dominate the 
transfer rate. Therefore, the increase of the efficiency 
and the pumping current with temperature is due to the 
increase of the number of diffusive trips of the shuttle. A 
temperature increase from 200 K to 400 K results in an 
increase of about a factor of two in the diffusion constant. 
It is expected that the proton current should increase at 
the same rate. However, our calculated ratio is about 
1.5 (see Fig. 10(a)). This is, probably, due to the fact 
that at high temperatures the shuttle has not enough 
time to be completely loaded with electrons and protons 
near the acceptor site A and the N-side of the membrane 
(and unloaded near the donor site D and the P-side of 
the membrane). A similar enhancement of the pumping 
current and the efficiency with temperature can be useful 
for photosynthetic microorganisms to compensate a leak- 
age of protons caused by the high-temperature increase 
of the membrane permeability |39j . The simple physi- 
cal features which come into play in our model arc also 
important for the creation of thermostable artificial pho- 
tosynthetic devices efficiently converting energy of light 
into electrical and chemical energy in a wide range of 
temperatures and light intensities. 
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FIG. 10: (Color online) (a) Proton current versus tempera- 
ture for different values of the light intensity I. (b) Pumping 
efficiency $ versus temperature. Here, the electrochemical 
gradient Afi = 220 meV (/Up = 110 meV, and fis = —110 
meV). 
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Effect of the electrochemical potential gradient 
on the proton current 



It follows from Eq. (3.3) that the difference, Afi = 
Mp — A*Ni between the electrochemical potentials of P- and 
N-proton reservoirs can be changed by changing the pH 
levels of the solutions inside and outside of the liposome. 
In doing so, one unit change in pH corresponds to ~ 
59 meV variation of the transmembrane proton gradient 
Afi (at standard conditions). To demonstrate the effect 
of the pH levels on the performance of the pump, in 
Fig. 11 we plot the dependencies of the proton current 
on the electrochemical potential ftp of the positive side 
of the membrane at three different values of the N-sidc 
potential: fi^ = -110; -140; and -200 meV. The proton 
current saturates when the P-side potential is sufficiently 
low, fip < 160 meV, and goes to zero at fip > 200 meV). 
At this condition, the potential of the P-side exceeds the 
energy, cq — 200 meV, of the proton on the shuttle: 
Hp > 6q, so that the proton cannot be translocated to 
the P-reservoir. On the other hand, the shuttle cannot 
be loaded with a proton at the N-side of the membrane 
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ber of diffusive trips of the shuttle. We also show that the 
pumping current increases linearly with the light inten- 
sity and saturates at the experimentally-observed limit, 
which is lower than the average intensity of solar light. 
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FIG. 11: (Color online) Proton pumping current versus elec- 
trochemical potential /ip of the positive side (P-reservoir) of 
the membrane for different values of the potential /ljn of the 
negative side (N-reservoir) for the light intensity / = 0.132 
mWcm~ 2 and temperature T = 298 K. 



if the electrochemical potential /zn is below the energy, 
6q — usq = —160 meV, of the proton on the shuttle 
populated with a single electron: p-^ < —160 meV. This 
is the reason why the last curve in Fig. 11 (taken at p,j$ = 
-200 meV) goes far below the other two curves (plotted 
for p^ > —160 meV). 



V. CONCLUSIONS 

We have analyzed a simple model for light-induced 
proton pumps in artificial photosynthetic systems. This 
model has five electron sites [four sites (D,B,C,A) for the 
triad molecule and one site for the shuttle (S)] and one 
proton-binding site on the shuttle (Q). The shuttle ex- 
hibits diffusive motion in the lipid bilayer, so that the 
electron and proton populations of the shuttle depend 
on the shuttle position. Based on the methods of quan- 
tum transport theory we have derived and solved numer- 
ically a system of master equations for electron and pro- 
ton state probabilities evolving in time together with the 
Langevin equation for the position of the shuttle. This al- 
lows us to calculate the proton current and the pumping 
efficiency of the system and determine their dependence 
on: the intensity of light, temperature and electrochem- 
ical potential gradient. 

For a reasonable set of parameters, closely related to 
the experimental setup, we demonstrate that this photo- 
synthetic device can translocate protons against an elec- 
trochemical gradient of the order of 220 meV with the ef- 
ficiency (photon-to-proton quantum yield) which exceeds 
55%. Our results explain experiments on artificial photo- 
synthetic reaction centers [10] . We predict that both the 
proton current and the pumping efficiency grow linearly 
with temperature due to the related increase of the num- 



APPENDIX A: BASIS SETS: 
ELECTRON-PROTON EIGENSTATES AND 
ENERGY EIGENVALUES 

The electron-proton system studied here with no leads 
can be characterized by the 20 basis states of the Hamil- 
tonian Hq 
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Here, |0) represents the vacuum state, when all electron 
and proton sites are empty. The state |1) = a^a^lO} 
corresponds to the case when one electron is located on 
the site D and one on the site B, and so on. The state 
1 11} = 44^q|0) indicates that, in addition to two elec- 
trons on the sites D and B, there is also a proton on the 
shuttle. The states |1) to |10) describe the shuttle with 
no protons, whereas the states |11) to 1 20) are related to 
the shuttle populated with a single proton. 

An arbitrary operator A of the combined electron- 
proton system can be expressed in terms of the basis 
Heisenberg matrices p m ,n — \m}(n\: 

m,n 

where m and n label the basis states: m,n = 1, . . . , 20. 
The diagonal operator is denoted as: p m = p m ,m- Thus 
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the electron population operators {ub, tib, n c, n A, «s} 
can be represented in the form 

TIB = Pl + P2 + P3 + Pi + Pll + P\1 + Pl3 + Pl7 

riB = Pi + Pa + P5 + P8 + Pll + Pi4 + Pl5 + Pl8 

nc = Pi + Pa + P6 + P9 + Pl2 + PlA + Pl6 + Pl9 

nA = P3 + P5 + P6 + PlO + Pl3 + Pl5 + Pl6 + P20 

ns = Pi + Ps + P9 + Pw + Pl7 + Pis + Pl9 + P20, 

(A2) 

and for the operator of the proton population of the shut- 
tle we obtain 

«Q = Pll + P12 + Pl3 + PlA + Pl5 + Pl6 + Pll + Pl8 + P19 + P20- 

(A3) 

Using the eigenfunctions (see Eq. (Al)), we can rewrite 
the Hamiltonian H in a simple diagonal form: 
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The terms 



J, 



(A5) 



describing the direct tunneling between 



all possible coupled sites i and i' are given by the expres- 
sions 



(A4) 



a B a D 

a B °c 
a^as 
a} D a s 



PA,2 + P5,3 + P8,7 + PlA.12 + Pl5,13 + Pl8,17 
P3,2 + P5.A + PlO, 9 + Pl3,12 + Pl5,14 + ^20,19 
Pl,2 + P5,6 + PH.9 + Pll, 12 + Pl5,16 + Pl8,19 
P3,7 + P5,8 + P6,9 + Pl3,17 + Pl5,18 + Pl6,19 
~Pl,8 — P2.9 — P3.10 — Pll, 18 — Pl2,19 — P13.20 

(A6) 



with the following energy spectrum: 



£1 


— Ey> + Eb 


; £11 


= £1 + £Q 


£2 


= Eb + Eq 


; £12 


= £2 + £q 


£3 


= E B + E A 


- u B a ; 


£13 = £3 + £Q 


£4 


= Eb + Eq 


- u B B ; 


£14 = £4 + 6q 


£5 


= Eb + E a 


- "da ; 


£15 = £5 + £Q 


£6 


= Eq + Ea 


- uba ; 


£16 = £6 + e Q 



It should be noted that the operator i/di r in Eq. (3) is 
non-diagonal. The proton operator bq can also be ex- 
pressed in a similar form 



bQ — J]] bQ^ nn p m , n . 



(A7) 
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